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0\ \ ABSTRACT 

'NT 

• We present synthetic Doppler maps of gaseous flows in binary IP Peg based on the 

' results of 3D gasdynamical simulations. Using of gasdynamical calculations alongside 

with the Doppler tomography technique permits us to identify the main features 
of the flow on the Doppler maps without solution of an ill-posed inverse problem. 
Comparison of synthetic tomograms with observations shows that in quiescence there 
^ ^, are two zones of high emission: a shock wave on the edge of the stream from Li 

I ' caused by the interaction of the gas of circumbinary envelope and the stream, and 

O I the dense region in a apoastron of quasi-elliptical accretion disk. A single arm of 

the spiral shock and the stream itself give a minor input to the total brightness. 



^ ' During outburst the accretion disk dominates, and the most emitive regions are the 

two arms of the spiral shock. 



spectrograms. Using this method one can ob- 
tain a map of luminosity in the 2D veloc- 



1 Introduction trograms for some emission line /(A, t) (or 

in other terms /(Vr, ip)) become widely used. 

„ J... 1 1 ri.- J. A method o£,Doppler tomography (Marsh & 

iraditional observations or bmary systems are ir>ool^\. , , , -i i 

. , ^ . , , i • 1 i Horne 1988c3 ) IS suited to analyze the trailed 

carried out using photometric and spectro- ' 

scopic methodics. The former gives the time 

dependence of brightness in a specific band . ^ , t ■ , . , -t r • 

rf,s 1,11,, • ,1 ,. 1 ity space trom the orbital variability ot emis- 

l{t) and the latter can give the time depen- rr^i , 

, „ , j_i r 1 1 1 sion lines intensity, i he Uoppler tomogram is 

dence ot wavelength ot some Doppler-shilted , . ^ , • i i 

\/,\ 1 • ■ 1 ii constructed as a conversion ot time resolved 

Ime Air). Given epnemeris Lp{t) is known, the ,. , r , , ,x ,. m ■ , ,i 

J J . rrj.\ J -L. i J (i-e. phase-tolded) line profiles into the map 

dependencies l(t) and X(t) can be converted ^ ,-,5^ Tr\ , r-^ , 

, rr^ 1 r i i ii i- ii On ( Kr , \/„ I plane, lo convert the distribution 

by virtue ot Doppler formula to the light curve \ • , t/tt ttn 

r/ \ J ii 1 J J r J- 1 HvRj'fi) into the Doppler map i(Vx,Vy) we 

7(031 and the phase dependency ot radial ve- , , , . r ,• , , 

T Tr , \ should use the expression tor radial velocity 

locity VR{ip). . ^. r 1 ^ T 

as a projection ot velocity vector on the line 

During the last ten years the observations of sight, i.e. Vr = -Va: cos{2Trip) + Vy sm{2Trip) 

of binary systems in the form of trailed spec- 
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(here we assume that 14 ~ 0, the minus sign 
before T4 is taken for consistency with the co- 
ordinate system), and solve an inverse prob- 
lem that is described by integral equatirin (see 
Appendix A of Marsh & Home 19880): 

^9{Vr + Vx cos(27r0) - Vy sm{2TT(j))) 
xdV^dVy , 

(1) 

where g{V) is a normalized local line pro- 
file (e.g., the Dirac (5- function) , and the lim- 
its of integration are from —oo to +oo. This 
inverse problem is ill-posed and the spe- 
cial regularization is necessary for its solving 
(e.g. Maximum Entropy Method, Narayaii_& 
Nityananda 19860, Marsh & Home 1988^3; 
Fourier Filtered Back-Projection, Robinson, 
Marsh & Smak 1993E31; Fast Maximum En- 
tropy Method— Spruit 1998c3l; etc., see also 
Frieden 1979ll3). 

As a result we obtain a map of distribution 
of specific line intensity in velocity space. This 
map is easier to interpret than original line 
profiles, moreover, the tomogram can show (or 
at least gives a hint to) some features of flow 
structure. In particular, the double-peaked 
line profiles corresponding to circular motion 
of the gas (e.g. in accretion disk) become a 
diffuse ring-shaped region in this map. Re- 
suming, we can say that components of binary 
system can be resolved in velocity space while 
they can not be spatially resolved through di- 
rect observations, so the Doppler tomography 
technique is a rather power tool for studying 
of binary systems. 

Unfortunately, the reconstruction of the 
spatial distribution of intensity on the basis 
of a Doppler map is an ambiguous problem 
since the points located far from each other 
may have equal radial velocities and deposit 
to the same pixel on the Doppler map. So the 
transformation /(Kc, Vy) I{x, y) is impossi- 
ble without some a priori assumptions on the 
velocity field. 

The situation changes drastically when one 
uses gasdynamical calculations alongside with 
Doppler tomography technique. In this 
case we don't need to cope with the in- 
verse problem since the task is solved di- 
rectly: p{x, y) & T{x, y) — » I{x, y) and 



I{x,y) & V,{x,y) & Vy{x,y) ^ I{V,,Vy) 
I{Vfi,(p). Difficulties can arise when convert- 
ing the spatial distributions of density and 
temperature p{x,y), T(x,y) into the distri- 
bution of luminosity of specific emission line 
I{x,y). For optically thick lines the forma- 
tion of line profile should be described by ra- 
diation transfer. equations (see, e.g.. Home 
& Marsh 1986c3), therefore to produce syn- 
thetic Doppler maps we assume that the mat- 
ter is optically thin and the intensity of con- 
sidered recomhipation line is as / ~ P^f~i^ 
(Ferland 1980EJ; Richards & Ratliff 1998E3). 

2 Binary system parame- 
ters 

The variable star IP Peg was discovered by 
Lipovetskij and Stepanyan (19810). It was 
found to be an eclipsing dwarf nova (or- 
bital period 3*^.79) with a deep eclipse and a 
humpjXm the light curve by Goransky et al. 
(19850). 

A binary system can be completely de- 
scribed by the following set of parameters: or- 
bital period P, masses of the components Mi 
and M2, separation A, and inclination angle 
i. Among all parameters only P is determined 
immediately from observations while param- 
eters Ml, M2, A, i (or in another statement 
M, q, A, i, where M = Mi + M2 - total 
mass, q — M2/M1 - mass ratio) should be 
calculated from observational data. These pa- 
rameters are connected by third Kepler's law 
n'^A^ = GM (here n ^ 2n/P - angular ve- 
locity of the orbital motion; G - gravitational 
constant), so we need only three additional re- 
lations. 

Several ways of determination of relation- 
ships between the parameters of the binary on 
the basis of observational data are known to 
exist: 

(i) using semi-amplitudes of radial velocity 

M2 q 
Ki = A - VL ■ sin i ■ — A-VL ■ sin i ■ , 

M 1 + g 



K2 = A-Vl ■ sin i ■ ^ A - VI ■ sin i ■ 

M 1 + q 

for mass-losing star and accretor, correspond- 
ingly, we can calculate q — K1/K2, and set a 
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relation between A and i (SI is considered to 
be already known) ; 

(ii) using the width of the white dwarf 
eclipse (for eclipsing binaries only) we can set 
a relation between q and i. 

Knowing observable entities Ki, K2, and 
A(j> permits us to determine all parameters of 
the binary system. Note, that alternate meth- 
ods of setting of relations between system's 
parameters also exist. For example, know- 
ing of rotational broadening of absorbtion lines 
Vrot = RRL{q, A)Q sin i, (where Rul ~ effec- 
tive radius of the Roche lobe) permits set an 
additional relation between A, q and i. Usage 
of extra relations can be exploited for checking 
of input parameters.]^ 

Semi-amplitudes of radial velocity for white 
dwarfjjjjere determined in. Wood & Crawford 
(19861^ ); Marsh (1988^); Marsh & Home 
(1990e3): Marsh (19880) obtained the value 
of Ki = 175 ±[15 km s^^ (but using Shaffer's 
method (1983H) he obtained a lower value 
of 164 km s-i). Wood & Crawford 

(1986E3) obtained the value of Ki = 141 
km s^^. Semi-amplitudes of radial velocity 
for mass-losing Jitar wcrc_jdctermined in Mar- 
tin et al. (19870, 198,90); Marsh (19^80); 
Beekmanetal (2000i): Marsh (19880) ob- 
tained the value 0L-K2 — 305 ± 15 km s~^, 
Martin et al. (19870, 19890) - the value of 
K2 =r?88 - 298 km s^i, and Beekman et al. 
(2000§) - the value ofi^2 = 331±5.8 km s~i. 
We adopt the values of semi-amplitudes of 
radial velotity according to Wood & Ccaw- 
ford (19860) and Marsh & Home (19900): 
Ki = 148 km s-\ K2 = 301 km s"^ 

After the values of Ki and K2 are adopted 
we can exploit the relation between system's 
parameters resulting from the value of width of 
the white dwarf eclipse. This value is known 
with sufficient accuracy: Wood & Crawford 
(19861 q) obtained Acj) = 0^0863, and Marsh 
(19880) obtained = 0^0858. The width 
of the white dwarf eclipse gives the relation 
between the orbit inclination i and the mass 
rati0j4,as (see, e.g., Home, Lanning &jComer 
19820; DhiUon, Marsh & Jones 19910): 



^Determination of system's parameters from width 
of 'hps spot' eclipsq— (see, e.g., Wood & Crawford 
IQSGH; Smak 1996Eal) is not considered here since 
we suggest an essentially different model of the flow 
structure in semidetached binaries. We also don't rely 
on the mass-radius relation since it is determined pri- 
marily for single stars. 
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Figure 1: Dependence between the orbit inclination 
i and the mass ratio q for the value of width of the 
white dwarf eclipse A0 = 0'\0863. A dashed line corre- 
sponds to an^pproximate formula jbii the Roche lobe 
radius (eq. (pi), see Pacynski 1971 123). A solid line 
corresponds to the exact formula for the Roche lobe 
sizes in Y and Z directions. Besides, we mrt values i 
and qp^opted by Wood & Crawfcyd (19860); Marsh 
(1988Ej); Beekman et al. (200otl) as weU as values 
adopted in our work: q = 0.49, i = 82°. 



A ■ tan(7rA0) • sin i 



RL 



A 



(2) 

Here ZjiL{q) and YRL{q) are the Roche lobe 
sizes in Z and Y directions (Z axis is perpen- 
dicular to the orbital plane and Y axis is di- 
rected along the orbital motion of mass-losing 
star). Usually ZfiL{q) and YfiLiq) functkms 
are taken approximately (Pacynski 197l'E3) 
as: 



ZRL{q)/A = YRL{q)/A = Rrl/A 



0.462 



1/3 



(3) 



Resulting dependency between i and q for the 
value A(j) = 0^0863 is shown in Fig. 1 (dashed 
line). However, usage of the approximate for- 
mula (H) may in this case lead to the loss of 
accuracy, so we have calculated dependencies 
ZRL{q) and YRLiq) exactly for 0.3 < q < 0.7 
and approximated it (error < 0.1%) as 

YRL{q)/A « -0.14858-(72+0.32031-g+0.18991 , 
ZRL{q)/A « -0.14466-g2-l-0.30559-q+0.18320. 
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Obtained relationship i(q) is shown in Fig. 1 
as weU (sohd hne). Besides, we put val- 
ues i jAt^d q adopted bii_,Wood & Crawford 
(19860)- Marsh (19880); and Beekman ct 
al. (2000§). 

Using the values for Ki, K2 mentioned 
above, as well as the more precise dependency 
i{q) we take the parameters for IP Peg as fol- 
lows: Ml = 1.02Mq, M2 = O.5M0, i = 82°. 
The distance between the inner Lagrangian 
point Li and accretor is D = 0.5 73 A = 
O.812i?0, the distance between system's center 
of mass and accretor is 0.329^1 — 0.466i?o. 

These parameters give the rotational broad- 
ening of the absorbtion lines of mass-losing 
star as Vrot — 143 km s^^. Observational 
measurements oLK-ot give Vrot — 146 km s"-'^ 
(Harlaftis 19990) and Vrot = fi^ km s^^ 
(Catalan, Smith & Jones 20010). Good 
agreement between calculated and observa- 
tional parameters in this additional relation 
bears witness on the correctness of system's 
parameters in use. 

We would like to stress that there were some 
factors noLincluded in our model. Martin^et 
al. (1987129) and Beekman et al. (2000i) 
point out a possible ellipticity of binary com- 
ponents' orbits (eccentrisity e — Q.QL — 0.075). 
Some autJbors (Wood et al. 1989H; Wolf et 
al. I993EI) argue that some observational pe- 
culiarities of IP Peg can be explained by the 
presence of the third body in the system. In 
our work we assume that the system contains 
two stellar components only, its orbits being 
circular. 



3 Model 

3.1 Gasdynamical model 

The full description of the 3D gasdynamical 
model can be found in Bisikalo et al. (2000O). 
Here we pay attention only to the main fea- 
tures of the model. To describe the gas flow 
in this binary system we used the 3D sys- 
tem of Euler equations for Cartesian coordi- 
nate system. To close the system of equa- 
tions, we used the equation of state of ideal 
gas with adiabatic index 7. To mimic the sys- 
tem with radiative losses, the value of adia- 
batic index has been accepted close to unit: 
7 = 1.01, that corresponds to the case close 
to the isothermal one (Sawada, Matsuda & 



Hachisu 198dc9 , Molteni, Belvedere. & Lan- 
zafame 19910; Bisikalo et al. 19950). 

The calculations were carried out in the non- 
inertial Cartesian coordinate system rotating 
with the binary system. The results of the 
calculations and the Doppler tomograms will 
be presented in coordinate system defined as 
follows: the origin of coordinates is located 
in the center of the accretor, X-axis is di- 
rected along the line connecting the centers of 
stars, from accretor to the mass-losing compo- 
nent, Y-axis is directed in the direction of or- 
bital movement of the donor-star, Z-axis is di- 
rected along the axis of rotation, so we obtain 
a right-hand coordinate system. Figure 2 (top 
panel) shows the coordinate system. In this 
figure we also put digits showing the phase an- 
gles of the observer in binary system, a Roche 
lobe with shadowed donor-star and the bal- 
listic trajectory of a particle moving from Li 
point to the accretor. The adopted coordi- 
nate system for Doppler maps is shown in the 
lower panel of Fig. 2. The transformation of 
the donor-star from spatial to velocity coor- 
dinate system is very simple as it is fixed in 
the corotating frame. Every point r fixed in 
the binary frame has a velocity x r in the 
corotation frame. This expression is linear in 
the perpendicular distance from the rotation 
axis, therefore the shape of the donor-star pro- 
jected on the orbital plane is preserved. Since 
the velocity of each point of the donor-star is 
perpendicular to the radius vector, all points 
of the donor-star are rotated by 90° counter- 
clockwise between the spatial and velocity co- 
ordinate diagrams (shadowed regions on the 
top and lower Panels of Fig. 2, see also Marsh 
& Horne 1988c3). On the velocity plane the 
accretor has coordinates (0,A'i). Figure 2 also 
depicts two concentric circles corresponding to 
different radii of the disk and their representa- 
tions on the velocity plane (for Keplerian ro- 
tation law). Inner circle has a larger velocity 
and forms the outer circle on the Doppler to- 
mogram. 

To obtain numerical solution of the system 
of equations we used the Roe-Osher TVD 
scherae of a high approximation order (Roe 
1986[Ea; Chakravarthy & Osher 19850) with 
Einfeldt modification (Einfeldt 19880). The 
computational domain was taken as a par- 
allepipedon [-D . . . D]x[D . . . D]x[0 . . . ^20] 
(due to the symmetry of the problem calcu- 
lations were conducted only in the top half- 
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Fi gure 2: Upper paneh The adopted coordinate system with phase angles of observer in a binary system. 
The asterisk is the accretor. Orbital rotation of the binary is counter-clockwise. The Roche lobe filling donor- 
star is shadowed. The critical Roche lobe and ballistic trajectory of a particle moving from Li are shown by a 
solid lines. Dashed and dash-dotted lines show concentric circles correspond to different radii of the disk. 

Lower panel: The adopted coordinate system in the velocity space. All designations are the same as in the 
upper panel. 
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space). A sphere with a radius of Vioo^ rep- 
resenting the accretor was cut out of the calcu- 
lation domain. The boundary conditions were 
taken as 'free outflow' on the accretor star and 
on the outer edges of computational domain. 
In the gridpoint corresponding to Li we in- 
jected the matter with parameters p = p{Li), 
T4 = c(Li), Vy = V:, = 0, where c(Li) is 
a gas speed of sound in Li point. Due to 
the scaling of the system of equations with re- 
spect to p (with simultaneous scaling of pres- 
sure P) we can accept an arbitrary value of 
p{Li) so we take it as p{Li) — l.Q The 
sound speed in Li was adopted as 5.5 km 
which corresponds to T{Li) = 3500 K. For 
the initial conditions we used rarefied gas with 
the following parameters po = 10~^ ■ p{Li), 
Po - 10-V(ii)c2(Li)/7, Vo = 0. 

3.2 Technique for construction 
of synthetic Doppler tomo- 
grams 

Preceding the consideration of the synthetic 
Doppler tomograms for binary system IP Peg 
it is necessary to stress the complexity of the 
analysis of these tomograms for eclinsiag sys- 
tems (see, e.g., Kaitchucket al. 1994IE3). The 
concept of Doppler tomography implies that 
any geometrical point deposits in some place 
on the velocity plane, this point being visible 
for any orbital phase. Indeed, from (|l]) we have 

I{Vr,^ + 7t)=I{-Vr,^), 

i.e. the constructions of synthetic Doppler 
tomograms is possible only for those sets of 
trailed spectrograms I(Vr, cj)) which are trans- 
formed symmetrically when one sees the bi- 
nary system from the "dark side" or, in other 
words, when there are no eclipses and occul- 
tations of emission regions. Clearly, eclips- 
ing systems violates that assumption. Usually, 
when dealing with observations, the eclipsed 
parts of trailed spectrograms arc naturally ex- 
cluded from input data for construction of 
Doppler tomograms. As a result, we obtain 
the Dopplcr map corresponding to the "trans- 
parent" case. The conversion of the results 

^When considering a system with known mass-loss 
rate, to determine the real values it is necessary to 
change the calculated values of density in accordance 
with the scale, defined by the ratio of the real value of 
the mass loss rate to the model one. 



of gasdynamical simulation into Doppler maps 
suggests using of full set of data that also cor- 
responds to the "transparent" case. 

The Doppler maps show the distribution of 
luminosity in the velocity space. Each point of 
the flow has a three-dimensional vector of ve- 
locity U = {Ux^ Uy, Uz) in observer's (inertial) 
frame. In the case when observer is located 
in the orbital plane of the binary the Doppler 
map's coordinate {Vx,Vy) will coincide with 
Ux and Uy. To define these coordinates for the 
case of inclined system we have to find a pro- 
jection of vector U on the plane constituted by 
vectors n and n x ft, where n is a direction 
from the observer to binary. 

The line emissivity in the velocity space can 
be written as: 



I{Vx,Vy)^ J J J I{x,y,z) 

O Uy 

'xS{Ux{x, y, z) sin i + Uz{x,y, z) coai — Vx) 
xd{Uy{x, y, z) sini -I- Uz{x, y, z) cos i — Vy) 
xdOdUxdUy , 

(4) 

where dO = dxdydz, i - inclination angle. 

As was mentioned above we adopt intensity 
as / ^ p'^T^/^ for the construction of synthetic 
Dopplcr tomograms. 

4 Results for quiescence of 
IP Peg 

Based on the model described in the Section 
4.1 we have conducted the 3D gasdynamical 
simulation of IP Peg in quiescence up to reach- 
ing of a steady-state solution. The morphol- 
ogy of gaseous fiows in considered binary sys- 
tem can be evaluated from Figs 3a. In Fig. 3a 
the distribution of density over the equatorial 
plane and velocity vectors are presented. In 
this Figure we also put a gasdynamical tra- 
jectory of a particle moving from Li to ac- 
cretor (a white line with circles) and a gasdy- 
namical trajectory passing through the shock 
wave along the stream edge (a black line with 
squares, see also Fig. 4). Analysis of the pre- 
sented results as well as our previous studies 
(Bisikalo et al. 19970, 1998bi) shows the 
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Fi gure 3cl! The distribution of density over the equatorial plane. Arrows are the velocity vectors in observer's 
frame. The asterisk is the accretor. The dashed-dotted line is Roche equipotential passing through L\. The 
white dotted line is the tidally induced spiral shock. Gasdynamical trajectory of a particle moving from L\ to 
the accretor is shown by a white line with circles. Another gas dynamical trajectory is shown by a black line 
with squares (see also Fig. 4). 



significant influence of tlie rarefied gas of cir- 
cumbinary envelope on the flow patterns in 
semidetaclied binaries. The gas of circumbi- 
nary envelope interacts with the stream of 
matter and deflects it. This leads, in par- 
ticular, to the shock-free (tangential) interac- 
tion between the stream and the outer edge 
of forming accretion disc, and, as the conse- 
quence, to the absence of 'hot spot' in the disc. 

At the same time it is seen, that the interac- 
tion of the gas of circumbinary envelope with 
the stream results in the formation of an ex- 
tended shock wave located along the stream 
edge ('hot line'). The 'hot line' model was 
confirmed by confronting j-with observations 
(see Bisikalo et al. 1998aD; Khruzina et al. 
200 lEJ). From Fig. 3a it is also seen the for- 



mation of tidally induced spiral shock (white 
dotted line in Fig. 3a). Appearance of the 
tidally induced two-armed spiral shock was nu- 
merically discovered in Saajada, Matsuda & 
Hachisu (1986a@, 1986b0\f,Sawada et al. 
(1987b3)j;-Spruit et al. (198713); Matsuda et 
al. (1990c3). Here we see only the one-armed 
spiral shock. In the place where the second 
arm should be the stream from L\ dominates 
and presumably prevents the formation of sec- 
ond arm of tidally induced spiral shock. 

An analysis of the flow structure out 
of equatorial plane shows that a part of 
the circumbinary envelope interacts with the 
(denser) gas stream and overflows it. This nat- 
urally leads to the formation, of 'halo'. Follow- 
ing to (Bisikalo et al. 2000O), one can define 
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Figure 3b: The distribution of p^T^/^ over the equatorial plane. A black dotted line is the shock wave along 
the edge of the stream ('hot line'). The main emission regions are marked by A, B, C, D. Other designations 
are the same as in Fig. 3a. 



'halo' as that matter which: i) encircles the ac- 
cretor being gravitationally captured; ii) does 
not belong to the accretion disc; iii) interacts 
with the stream (collides with it and/or over- 
flows it); iv) after the interaction either be- 
comes a part of the accretion disc or leaves 
the system. 

Figure 3b depicts the distribution of p^T^f^ 
over the equatorial plane. Similar to Fig. 3a 
spiral shock is shown by white dotted line. Be- 
sides, shock wave along the edge of the stream 
is shown by black dotted line. The distribu- 
tion shown in Fig. 3b represents the intensity 
of recombination line, so the analysis of this 
distribution can determine the most luminous 
region of the flow. It is seen, that the main 
emission comes from four region designated by 
markers A, B, C, D. 



• Marker A designates the shock wave 
along the edge of the stream ('hot line') 
resulting from the gasdynamical interac- 
tion of the gas of circumbinary envelope 
with the stream. 

• Marker B designates the stream from Li 

or, more exactly, the most luminous part 
of the stream where the density is still 
large enough and the temperature already 
increases due to dissipation. 

• Marker C designates a region near the 
apoastron of the accretion disk. The anal- 
ysis of the presented results shows that 
the disk has a quasi-elliptical form, there- 
fore approaching the apoastron the mat- 
ter is retarded and the dense region is 
formed. 
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Figure 4: Synthetic Doppler map for I ~ p^T^I"^. The secondary Roche lobe (a bold black line) and the 
accretor (an asterisk) are also shown. The white line with circles and black line with squares show gas dynamical 
trajectories in the velocity coordinates (see Fig. 3a). The main emission regions are marked by A, B, C, D as 
in Fig. 3b. 



• Marker D designates a dense post-shock 
region attached to the spiral shock. 

The synthetic Doppler map based on the re- 
sults of 3D gasdynamical simulations is pre- 
sented in Fig. 4. Earlier we have analyzed the 
features of flow in the equatorial plane (see 
Figs 3a, 3b), but for the sake of comparison 
with observations synthetic Doppler maps will 
be presented for integrated over z-coordinate 
intensity in accordance to the equation (jj). 
Two gasdynamical trajectories (the same as in 
Fig. 3a but in velocity coordinates) are shown 
in Fig. 4. 

Shock wave resulting from the gasdynami- 
cal interaction of the gas of circumbinary en- 
velope with the stream is located along the 
edge of the stream. Three last points (marked 



by larger symbols) of curves with circles and 
squares in Fig. 3a are the examples of two flow- 
lines passing through the shock. Location of 
these parts of trajectories on the Doppler map 
corresponds to region A left to the donor-star. 
This region of Doppler map contains also a 
spiral arm beginning approximately from the 
center of mass-losing star and located above 
the vicinity of L\. Our analysis shows that 
the appearance of this region results from the 
overflowing the stream by the gas of circumbi- 
nary envelope but not from the shock wave. 
It is seen from Fig. 4 that the stream from 
L\ (the beginning of white line whit circles) 
transforms into spiral arm B in III quadrantn 

^The quadrants of coordinates plane are counted 
as follows: I quadrant is upper right (corresponding to 
Va; > 0, Vj/ > 0), other quadrants are counted counter- 
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Figure 5: Doppler maps of H-y and lines for IP Peg in quiescence (Wolf et al. 19980). This figure is 
reproduced under the kind permission by S.Wolf and A.Bobinger. 



of the Doppler map. The region of increased 
density near apoastron of the disk - region C 
is seen in Fig. 4 as more luminous zone on the 
border of I and IV quadrants. Tidally induced 
spiral shock (or, more exactly, the dense post- 
shock zone, dotted line in Figs 3a, 3b) forms 
a bright arm in I and II quadrants of Doppler 
map. 

Resuming these results for Doppler map of 
IP Peg in quiescence we can conclude that 
there are four elements of the flow structure 
which deposit in the total luminosity: 'hot 
line', the most luminous part of the stream 
where, the dense region near the apoastron of 
the disk, and the dense post-shock region at- 
tached to the spiral shock. The income of each 
element obviously can vary depending on pe- 
culiarities of considered binary system. It is 
also obvious that based on the model compu- 
tations we can't estimate what elements will 
dominate. Nevertheless, comparison of syn- 
thetic Doppler maps and observed ones per- 
mits both to catch the dominating element and 
to correct/refine the computational model. 



Observational Doppler tomograms for IP 
Peg in-quiescence were built in Matah & Horne 
(1990El)i_Harlaftis et al. (1994113 )■ Wolf et 
al. (1998[|lL Bobinger et al. (1999i); Bob- 
inger (2000ltJ). Figure 5 represents a typical 
Doppler map for and from Wolf et al. 
(1998H). The characteristic features of these 
tomograms are the bright spot in the region 
A as well as the zone of moderate brightness 
in the region C. The comparison of the obser- 
vational tomogram from Fig. 5 and synthetic 
one from Fig. 4 reveals that it is the 'hot line' 
and the dense zone near the disk's apoastron 
which mainly deposit into the total luminos- 
ity. Signatures of the spiral shock are not seen 
and this implies either its absence or weak- 
ness. Note also that the observational tomo- 
gram shows rather small input from the stream 
from Li into the total luminosity. 



clockwise. 
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Fi gure 6(1: The distribution of density over the equatorial plane for outburst. Arrows are the velocity vectors 
in observer's frame. The asterisk is the accretor. The dashed-dotted line is Roche equipotential passing through 
Li . White dotted lines are the arms of tidally induced spiral shock. 



5 Results for outburst of 
IP Peg 

Observations show (see, e.g., recent— reviews 
by Marsh 2000^ and Steeghs 200oll3) that 
during outburst the accretion disk dominates 
hence the stream from L\ plays less impor- 
tant role. Our today's knowledge of the na- 
ture of the outburst as well as its parame- 
ters has an approximate and qualitative char- 
acter so it is hard to simulate the outburst 
correctly. To mimic the flow structure dur- 
ing outburst on the qualitative level we cal- 
culated the structure of gaseous flows up to 
reaching a quasi-steady-state solution and put 
the rate of mass transfer equal to zero (i.e. ter- 
minated the mass transfe|t) as it was suggested 
in Bisikalo et al. (2OOla0, 2001bi). We un- 



derstand that this model doesn't reflect all pe- 
culiarities of outburst and expanding accretion 
disk but we hope that it truly correlates the 
influence the disk and the stream on the qual- 
itative level. Our simulations of residual ac- 
cretion disk show that at time 0.3 QAPorh 
after mass transfer termination the flow struc- 
ture is changed significantly. The stream from 
L\ vanishes and doesn't dominate anymore, 
and the shape of accretion disk changes from 
quasi-elliptical to circular. The second arm 
of tidally induced spiral shock is formed while 
earlier (before the termination of mass trans- 
fer) it was suppressed by the stream from L\. 
It is seen that obtained flow structure has all 
basic features observed in outburst of IP Peg. 
This gives a hope that we can refine/reveal the 
new features of IP Peg in outburst by virtue of 
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Figure 6b: The distribution of p^T^/-^ over the equatorial plane for outburst. The main emission regions are 
marked by A and B. Other designations are the same as in Fig. 6a. 



analysis of synthetic Doppler tomograms con- 
structed for this gasdynamical solution. 

Figure 6a depicts the distribution of den- 
sity and velocity vectors and Fig. 6b depicts 
the distribution of p^T^^^ over the equatorial 
plane. It is seen, that the main emission comes 
from two arms of the spiral shock designated 
by markers A and B. The synthetic Doppler 
map based on the results of 3D gasdynamical 
simulations for IP Peg in outburst is presented 
in Fig. 7. Our analysis shows that bright arms 
in I and III quadrants are due to emission of 
dense post-shock zones attached to the arms 
of spiral shock. 

Observational Doppler tomograms for IP 
Peg during the Jiutburst were built in Maxsh 
& Home (19900); Steeghs et al. LL9960); 
Steeghs, Harlaftifj_f: Home (1997H); Har- 
laftis et al. (1999E3); Morales- Rueda, Marsh 



& BiUington (20000). A typical example 
of such tomogram (Morales-Rueda, Marsh & 
BiUington 2000^) is given in Fig. 8. The 
characteristic features of these tomograms are 
two bright arms in I and HI quadrants. The 
comparison of the observational tomogram 
from Fig. 8 and synthetic one from Fig. 7 re- 
veals that these arms results from dense post- 
shock zones attached to the arms of spiral 
shock (zones A and B). 



6 Conclusions 

Using of the gasdynamical calculations along- 
side with Doppler tomography technique per- 
mits us to identify main features of the flow 
on the Doppler maps without solution of the 
ill-posed inverse problem. The comparison 
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Fi gurc 7; Synthetic Doppler map for the outburst. The secondary Rx)che lobe (a bold black line) and the 
accretor (an asterisk) are also shown. The main emission regions are marked by A and B as in Fig. 3b. 



of synthetic Doppler maps and observed ones 
permits to correct/refine the computational 
model and to interpret the observational data. 
In this work we have presented the synthetic 
Doppler maps of gaseous flows in binary IP 
Peg based on the results of 3D gasdynami- 
cal simulations. The zones of flow structure 
responsible for the most emitive regions of 
Doppler map were identified and it was found 
that they are different for quiescence and out- 
burst. Our analysis for quiescence has shown 
that it is the shock wave along the stream edge 
- 'hot line' and the dense zone near the disk's 
apoastron which mainly deposit into the total 
luminosity. The input from the stream from 
Li and the spiral shock into the total lumi- 
nosity is small. During the outburst the role 
of the stream is unimportant and the accre- 
tion disk with two-armed spiral shock domi- 



nates. The comparison of the observational 
tomogram and synthetic one reveals that the 
bright arms in the Doppler map result from 
dense post-shock zones attached to the arms 
of the spiral shock. 
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